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Abstract 

The three-dimensional boundary layer on a swept wing can support different types of 
hydrodynamic instability. Here attention is focused on the so-called ‘spanwise contamina- 
tion’ problem which occurs when the attachment line boundary layer on the leading edge 
becomes unstable to Tollmien-Schlichting waves. In order to gain insight into the inter- 
actions which are important in that problem a simplified basic state is considered. This 
simplified flow corresponds to the swept attachment line boundary layer on an infinite flat 
plate. The basic flow here is an exact solution of the Navier Stokes equations and its sta- 
bility to two-dimensional waves propagating along the attachment line can be considered 
exactly at finite Reynolds number. This has been done in the linear and weakly nonlinear 
regimes by Hall, Malik and Poll (1984) and Hall and Malik (1986). Here the corresponding 
problem is studied for oblique waves and their interaction with two-dimensional waves is 
investigated. In fact oblique modes cannot be described exactly at finite Reynolds number 
so it is necessary to make a high Reynolds number approximation and use triple deck 
theory. It is shown that there are two types of oblique wave which, if excited, cause the 
destabilization of the two-dimensional mode and the breakdown of the disturbed flow at 
a finite distance from the leading edge. Firstly a low frequency mode closely related to 
the viscous stationary crossflow mode discussed by Hall (1986) and MacKerrell (1987) is a 
possible cause of breakdown. Secondly a class of oblique wave with frequency comparable 
with that of the two-dimensional mode is another cause of breakdown. It is shown that 
the relative importance of the modes depends on the distance from the attachment line. 


* This research was supported by the National Aeronautics and Space Administration under 
NASA Contract No. NAS1-18107 while the authors were in residence at the Institute for Computer 
Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, 
VA 23665. 


1. Introduction 

Recent interest in the development of laminar flow wings has generated much interest 
in the different modes of instability which the three- dimensional boundary layer on a swept 
wing can support. There are essentially three different types of instability which such flows 
can support. Firstly these flows are susceptible to Tollmien-Schlichting instability waves 
induced by viscosity. These axe important everywhere in the flow. Secondly the crossflow 
instability mechanism discussed by Gregory, Stuart and Walker (1955) is operational in 
three-dimensional boundary layers. This mechanism is inviscid in origin and is important 
sufficiently far from the leading edge of the wing where the spanwise velocity component 
of the flow is not negligible. Finally, if there are regions of concave curvature on the wing, 
the centrifugal instability mechanism will lead to counter-rotating Gortler vortices. 

Clearly in some boundary layers all three mechanism will be operational and their 
interaction might induce the premature transition of the flow. However, it is clear that if 
the boundary layer does not remain laminar in the attachment line region then there is no 
point in trying to control the other mechanisms further away from the leading edge. 

Experimental investigations of the attachment line instability problem have been made 
by Pfenninger and Bacon (1969), Gaster (1967) and Poll (1979). It was found that, if the 
input disturbances are small, instability waves propagate along the attachment line when 
the flow Reynolds number exceeds a certain critical value. Pfenninger and Bacon (1969) 
found that if the input disturbances were sufficiently large then there was some evidence of 
a subcritical response. The measurements made in these experiments were all close to the 
leading edge and showed no evidence of oblique Tollmien-Schlichting waves being present 
in the flow. 

Here we shall first discuss the linear instability of the attachment line boundary layer 
on an infinite flat wall. We describe the instability using the approach of Smith (1979a, b) 
and therefore use triple deck theory. We find that in addition to the two-dimensional 
mode of Hall, Malik and Poll (1984) (hereafter referred to as HMP) there is possible, at 
any chordwise location, a family of unstable oblique waves. Each oblique mode is unstable 
only up to a finite value of X, the chordwise variable. This instability problem is discussed 
in §2. 

We then consider the nonlinear interaction of a two-dimensional mode with an unstable 
oblique mode. This is done using the formulation of Hall and Smith (1984) for wave 
interactions in Blasius flow and leads to a pair of coupled amplitude equations which 
determine the evolution of the modes. These equations are found in §3 and in §4 we discuss 
the possible solution of these equations for different input waves. We show that there is a 
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particular type of interaction which leads to both disturbance amplitudes becoming infinite 
at a finite chordwise location. 

In §5 we briefly describe the numerical work needed to calculate the coefficients in 
the amplitude equations. We then discuss the nature of the solutions of the amplitude 
equations appropriate to the constants we calculate and the implications for the stability 
of the flow of interest. 


2. Formulation of the stability problem 

Consider the flow of a viscous incompressible fluid of kinematic viscosity v adjacent 
to the flat plate defined by y = 0 with respect to Cartesian co-ordinates ( x,y,z ). The 
velocity field (u,u,u>) corresponding to (x, y, z) satisfies the conditions 


u = i> = u> = 0, y = 0, 

X 

u->Uooj, w->Uoo, y — ► oo. 

The Reynolds number R is defined by 


R = 


Uooi 




V 

and dimensionless variables ( X , Y, Z ), (U, V, W) are defined by 


(2.1a, b ) 


( 2 . 2 ) 


(X,Y,Z) = ( 

0 u,v,w) = U„(U,V,W ). 


(2.3a, 6) 


If the pressure is scaled on pU^ and time, T, on Uoo/£, then the Navier Stokes equations 
become 

U T + (U* V)C7 = - VP + 4v 2 P. (2.4) 

R 

In order to satisfy the no-slip condition the basic state has a boundary layer of thickness 
£R~i near Y = 0. If the boundary layer variable rj is defined by 


rj = R^ Y, 


the basic state can be expressed in the form 

U = (XuM , R-ivMMvm + O(R-i)), (2.5) 

where 

u + v = 0, 
v'" + v' 2 — vv" — 1 = 0, 
w" — vw = 0, 


(2.6a, b, c ) 
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with 


u(0) = u'(0) = u>(0) = 0, 
v'(oo ) = — 1 , w( oo) = 1 . 


(2.7a, b) 


It is known, Smith (1979a, b), that lower branch Tollmien-Schlichting instabilities of bound- 
ary layers are governed by triple deck theory. Here the interest is in oblique Tollmien- 
Schlichting waves proportional to E where 


E = exp i 




cO cl 


( 2 . 8 ) 


with e = R s << 1. The slowly varying wavenumber a is then expanded as 


a = cto + eoq H . 


(2.9) 


In the main part of the boundary layer the basic state is perturbed by writing 

U = (Xu, e 4 v,w)( 1 + 0(e 4 )) + 6(eU 0 ,e 2 V 0 , eW 0 )E + • • • , (2.10) 


where Uq , Vo and Wo depend only on X and rj whilst 8 is assumed small compared to any 
power of e. The corresponding pressure perturbation is e 2 P 0 E where Po is a function of X 
only. The equations to determine Uq , Vo in the main deck are: 


iaoUo + Vo,) + */?Wo = 0, 

iaoXuUo + XVqv! + i{3wUo = 0, (2.11a, b, c ) 

iaoXuWo + Vqw' + i(3wWo = 0 , 

and the appropriate solution of this system is 

U 0 = a(X)Xu ' , Wo = a(X)u;', 

(2.12a, b ) 

V 0 = -ia(X)[a 0 Xu + (3w], 


where a(X ) is an amplitude function to be determined. An investigation of the disturbed 
flow in the upper layer shows that all disturbance quantities decay exponentially there and 
that matching with the main solution requires 

Poyf^T7 2 = + /?] 2 . (2.13) 


It should be noted from (2.12) that when // — * 0 


U 0 — a(X)AAT , W 0 ^ a(X)n , 


(2.14) 
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where from numerical calculations it is found that A = 1.236, [x = 0.570. Finally in the 
lower deck the pressure perturbation is still € 2 PqE whilst the total velocity field can be 


written : 

U = e£AX + • • • Scu 0 E 
V = -e 6 ( 2 \X/2 + --- +6e 3 v 0 E + 

W = + • • • + 8ew 0 E + • • • . 

Here the lower deck variable C, is defined by 


(2.15a, b, c ) 


( = Ye~ 5 . 


The equations to determine (u 0 , vo> wo) the l° wer deck can be written as 


C V ° = 0 , 

Wo 


(2.16) 


where the matrix operator 


— iaoXA( — iftixC + —AX 


— n jgz — iaoXA( — i0nC + — i/3 


The solution of the wall layer problem is then written as 


OcqUq + (3w 0 


where 


- /3wo = b(X ) / Ai(s)ds , 

Jio 

Avo = ~i [tto^o + /3wo\ds, 

Jt 0 

A = {iAX<* 0 + ifii*}*, 

e = Ac+e 0 , 

iQ, 

= ' 


(2.17) 


(2.18) 


(2.19) 


( 2 . 20 ) 


The pressure Pq can be related to the displacement function b using the X and Z momen- 
tum equations to give 

i{a 2 0 +p 2 }P 0 = A 2 Ai'(h)b. (2.21) 
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The lower and main deck solutions are then found to match if 

Ax' o)(<* 0 AX + /fy) 1 = iiy/al + P X o(<*oX + 0)\ (2.22) 


where _ 

pOO 

Xo = / At’(a)ia. 

Ao 

The eigenrelation (2.22) determines the complex wavenumber ao for given 0 and f l. For 
neutral stability it is known that 

£o ^ —2.298*5, ~ 1.001*4 . (2.23a, b) 

Xo 

Thus the neutral values of ao and 0 are related by 

1.001(a 0 AX + /3fi)s ~ ypi + /? 2 (a 0 X + 0f. (2.24) 

The modes corresponding to HMP have ao = 0 in which case 0 satisfies 

0* ~ 1.001 /J. 

At any given value of X there are in addition neutral three-dimensional modes with ao ^ 0. 
The above analysis fails if a 0 X + 0 or aoAX + 0(x become negative anywhere in the (a o ,0) 
plane so only eigenvalues of (2.24) above the lines ao AX + 0fj, = 0, and aoX + 0 = 0 are 
acceptable. In Figure la ao is shown as a function of 0 for X = 0.1,1., 10., 20., 30.. The 
solutions in the second quadrant asymptote to the line aoX + /? = 0 as /3 — * — oo whilst 
those in the fourth quadrant asymptote to the line ao AX + 0[i = 0 as 0 — » 0. Figure lb 
shows the neutral values of as a function of 0 for X = 0.1, 1., 10., 20., 30.. 

Finally in this section we notice that the two-dimensional mode of (2.22) which of 
course corresponds to ao = 0 is neutrally stable at all values of X. The three-dimensional 
modes however axe initially unstable on the attachment line X = 0 and become stable 
beyond a critical value of X. Experimentally it appears that if the level of disturbances 
present in the flow is sufficiently small then it is the two-dimensional mode which is ob- 
served. In the next section we investigate the possibility that the two-dimensional mode 
might be destabilized by oblique modes which grow in the X direction. In Figure 2 we 
have shown typical growth rate curves for the three-dimensional modes. 

In fact the small 0 solutions are related to the stationary modes of instability of the 
three-dimensional boundary layer discussed by Hall (1986) and MacKerrell (1987). These 
modes orient themselves such that the shear stress of the ‘effective’ velocity profile is zero; 
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the lower deck structure is then described by parabolic cylinder functions rather than Airy 
functions. Thus when Oo tends to zero we find from (2.24) that 

/i/3 = — c*o^AT + 0(o;o)® > 

6 

so that the neutral frequency tends to zero like o 0 5 . Some discussion of the time dependent 
version of the stationary modes discussed by Hall and MacKerrell has recently been given 
by Bassom (1987). 

3. Weakly nonlinear theory 

Suppose that the three-dimensional mode with (a,/3,0) = ( 0 : 2 , P 21 CI 2 ) is neutrally 
stable at X = X n . We consider the interaction of this mode with the two-dimensional 
disturbance which propagates along the attachment line. We know from the work of Hall 
and Smith (1984) that in the absence of the two-dimensional mode the three-dimensional 
mode will evolve in a nonlinear, nonparallel manner in an ef neighbourhood of X n . We 
therefore define X by 

* = (3.1) 

Later we can derive the ‘quasi-parallel’ evolution equations for (X — X n ) > 0(e^ ) by taking 
the limit X — ► 00 . In order that the two-dimensional mode in this neighbourhood should 
be of finite size we suppose that, with Cl = Cli, the neutral frequency for a two-dimensional 
wave, the spanwise wavenumber pi is expanded as 

Pi — P 10 + €/?n + e 2 P + • • • , (3-2) 

where Pio,Pu are the first two terms in the expansion of the neutral spanwise wavenumber. 

It is now convenient to represent the ‘fast’ dependence of the Tollmien-Schlichting 
waves in the X direction by multiple scales rather than the WKB formulation of §2. We 
therefore write 

X * = - X - ~ Xn ^ (3.3) 

r.rPlZ 

Ei =exp[*{-3 —}], 

ITT r-r ta* /?2 Z CI 2 T 

E 2 = exp[*{a 2 A: + — —}], 

e 5 e 

1 

\ 


Next we define 
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where fti , O2 and /? 2 expand as 

fti = ft 10 + eft 11 + 0(e 2 ), 

ft2 = ft20 d* eft2i -f" 0(e 2 ), (3.4a, 6, c) 

/?2 = P 20 + e/?2i + ^( e2 )- 

Here ftio,ftn etc. are the neutral values appropriate to the location X = X n whilst j 3 \ is 
as given by (3.2). For the three- dimensional mode we further expand 

a 2 = ®2o ea2i + 0 (e 2 ), (3.4 d) 

where 020, <*21 are the first two terms in the expansion of the neutral value of a 2 at X = X n . 
In the lower deck we write the velocity in the form 

U = e{\X + • • • +eU, 

V = -e 6 C 2 AX/2 H +e 3 V, (3.5) 

W = eC,n + eW, 

and then expand the disturbance velocity field (U, V, W) together with the corresponding 
pressure perturbation as 

(U, V, W, P) = e * & + e* S 2 + e* S3 + el S* + • • • . (3.6) 

Here the term S_i corresponds to the fundamental modes proportional to E\ and E 2 - 
The second order term S 2 corresponds to first harmonic and mean flow correction terms 
generated by the interaction of the TS waves. The third order term S3 again contains the 
fundamentals generated because the correction terms in (3.4) are 0(e). Finally the fourth 
order term S ' A contains fundamental and other terms driven by the interaction of S3 and 

Clearly the function S3 satisfies the linearized problem of §2 so we write 

5j = AS^Ex + BS 12 E 2 + C.C. (3.7) 

where C.C. denotes ‘complex conjugate’ and A, B axe functions of X to be found at higher 
order. The functions S_ij are defined by S_ij = ( U t j , Vij , W{j , P,y ) for i,j > 1 and S n , S ]2 
satisfy (2.16) with (ao,/ 3 , ft) replaced by (0, /?io, ftio) and (020^20^20) respectively. Thus 
for example we can show that 
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<* 20^12 + 02 oW \ 2 — / Ai(s)ds, 

J $20 

6o = --^» £2 = A 2 C + £ 2 o, A 2 = {iAJ„a 2 o + */? 2 om} 3 ' 

At next order we find that the first harmonics and mean flow corrections can be 


written as 


s 2 ={A 2 S 21 E 2 + B 2 S 22 E 2 + ABS 23 E 1 E 2 + ABS^EiE- 1 } + C.C. 
+ | A \ 2 S 25 + | B | 2 S 26 , 


(3.8) 


where B denotes the complex conjugate of B. We find that S 21 , S _ 22 satisfy the differential 
system 


£(2a n o, 2/? n 0} 2f2 n o) 


( V2n \ ( ia n 0 U 2 n + V ln ^ + i/3 n 0 W ln U l7 

V2n 

w 2n 
\P2n I 


d< 

0 


,ia n 0 U ln W ln + V ln ^ + i0 nO W? n 


(3.9) 


for n = 1,2 and aio = 0. These equations must be solved subject to 


u 2n = v 2n = w 2n = 0 , C = 0 . 


(3.10) 


The functions S 23 and S 24 satisfy similar equations but with (2a n o, 2/? n o, 2fl n o) re- 
placed by (c*io ± CH 20 , 0io ± 020, ftio ± ^ 20 ) for n = 3,4 and the right-hand-side of (3.9) 
replaced respectively by 

ioc 2 oU \\Ui2 + V \ 2 + Vu dJ ^,\ 2 + i[0ioW 12 Un + 0 2 oWnUi 2 ] 

0 

ia 20 UuW 12 + V 12 ^ + Vn^ + i[0 w + 0 2 o]W u W 12 


and 


-^20^11^12 + V12 -3^ + Vn + i[/?io W\iU\i — ^20^11^712] 

0 

-i a20 UiiWi 2 +Vi 2 ^ + V 1 i^ + i(0 lo -02o)W u Wi2 

The mean flow corrections for 5 25 and S _ 26 have V 25 = V 2 e — 0 whilst for n = 1,2 


d f^*2n-f*4 jr dUin dUin /... __ jj . 

^q2 — 'In A 'In T %0nO\ '' InBln ^^lnUln), 

d W2n+4 /fj jjf TT j~jj r ^ , T/ dWjn ( ~ r dW\ n 

— l&nOyEln Win b^lnWi n ) + Vj n , A *b Tin 


(3.11a, 6) 


d c 2 


dC 


dC ’ 
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which must be solved subject to 


U 2 n-\-4 Hdn+4 — 0, C — 0, 

dU 2n+ 4 <»Wi „ , (3.12a, 6) 

“^ = ~^ = °’ C = 0 °' 

In the main and upper decks the mean flow corrections and the first harmonic func- 
tions are not forced by the fundamentals and therefore satisfy similar equations to those 
discussed in §2; the matching of the main deck solutions for the first harmonics with the 
lower deck solution produces boundary conditions at £ = oo for 5' ?1 ,5' ?? ,5 , 2 ? ,5 , ?4 . 

At next order in the lower deck problem we obtain only fundamental terms driven 
by the variation of the mean state. In fact the solution of this linear problem when 
matched with the main deck solution determines the 0(e) terms in the expansion of the 
neutral wavenumber and frequencies. Since the solution at this order has no effect on the 
amplitude equations for A and B we give no details of it here. 

The interaction of the fundamental term S_ x with the mean flow correction and first 
harmonic term generates fundamental terms in S A . In addition further fundamental terms 
are produced by the evolution of the amplitude functions A and B and the basic state on 
the X length scale. If we write in the form 


*£4 — *£41 i 5 i + >£42^2 + C.C. + • • • , 


where • • • represents other terms forced by the interactions, then after some manipulation 
we find that SUi and •£42 satisfy 


m /?io, ^10) 


/U 41 

V 41 

W 41 


VP 41 


20, # 20 ,^ 20 ) 


U 42 \ 

v 42 

W 42 

V P 42 / 


[A | A | 2 + A | B | 2 ^2+ifiPCUuA \ 

+\XVnA + \X n (Uujj£ + Pn ' 

-Unj^-iWuPA 

[A | A | 2 $ 3 + A | B l 2 $ 4 + ifoCWnA 
+\CX n W 11 j$ + 0 P 11 A] / 


| 2 $s + P | A | 2 $6 + AX„Cg| U 12 \ 
n X\CBU,o + \XBV10 -l-P,o^l * 


\[B | B | 2 $ 7 + B | A | 2 $ 8 + AC X n W 12 ^ 

\ +\(lXWi 2 ia 2 QB] J 


Here the functions $ 1 , $ 2 , $3 and $ 4 are defined by 


$1 = iaio(U 25 Un+U 21 U 11 )+V 11 U^+V 21 U[ 1 +VnU^+i 0 lo (W 25 U 11 -W 21 Un+ 2 W n U 21 ), 
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$2 =*«lo(^12^24 + U 2 e,U\\ + U 23 U 12 ) + V\2U' 2i + V 23 U 12 + ^24^12 

+ ^12^23 + ^11^26 + *(£l0 + ^20)^12^23 + i{Pl 0 — /?2o)Wl2^24 
+ ^ 20 ^ 24^12 — 2 /^ 20 ^ 23^12 + 2 /^ 10 ^ 26 ^ 11 , 

$3 =ia 10 (U 25 Wu + 2tf n W 21 - tf 2 iWn) + V„W' 5 + V 2 iWu 

+ Vn W 2 \ + 2^10(^25^11 + ^11^21), 

$4 =2(o!lo + <^20)^12^23 + 2(0110 — 020)^12^24 — 2020^23^12 

+ 2020^24^12 + 2010^26^11 + + V23W12 + ^24^12 (3.13a, b , C, d ) 

+ V 12 W^ + VuW ^ 6 + ip 10 (W 12 W 24 + W 26 W n + W23W12), 

with 010 = 0. The corresponding expressions for $5, $6, $7, and $8 are obtained from 
(3.13) with (010, /?io,^io) and (020,^20,^20) interchanged and then 010 set equal to zero 
together with 

{A, Bsuffixes 11, 12, 21, 22, 23, 24, 25, 26} 


replaced by 

{B, ^suffixes 12, 11, 22, 21, 23, 24, 26, 25}. 

Finally the terms with suffix 24 are replaced by their complex conjugate. The disturbance 
velocities (U41, V41, W41) and (U 42 ,V 42 ,W 42 ) must of course vanish at £ = 0 and the 
functions ^415^2 must match with the corresponding functions in the main deck. The 
latter matching conditions completely specify inhomogeneous differential systems for £41 
and §42- Since the homogeneous form of these systems have a solution it follows that we 
must apply solvability conditions to the systems for ^41>^42- In order to write down these 
conditions we must introduce the differential systems adjoint to those which determine the 
fundamentals. We first note that if we define 


F = a Q u 0 + Pwo, 


in (2.18) then the eigenvalue problem which leads to (2.22) can be written as 

F"' - f[a 0 AX + PfACF' + itlF' = 0, 

F( 0) = F'(oo) = 0, (3.14) 

F(oo) = CF"{ 0), 


-i[q 0 AX + Pfj] 

(<4 + P)k( a „x + ay' 


where 
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The system adjoint to (3.14) is 

p' = 0 , 

r" — i[aoAX + — iXlr = p, 

g( 0 ) = r(oo) = 0 , r( 0 ) = Cp(oo). 


(3.15) 


It is easily seen that we can take p = 1 above and then solving the equation for r using 
variation of parameters and the boundary conditions give (2.22) again. It then follows that 
if f2 = Q,(ao,/3) is an eigenvalue of (3.14) then the system 

G'" - i(a 0 X A + + ittG' = R, 

G( 0 ) = G'( oo) = 0 , G(oo) - CG"( 0 ) = 7 , 

will have a solution if 

poo 

/ rRd C = 7 . 

Jo 

Thus it follows that the differential systems for S A1 ,S A2 will have a solution if 

dA 


dX 

dB_ 

dX 


= X^A-aU | A | 2 -hA | B |' 
= A 2 XB - a 2 B | B | 2 -b 2 B \ A 


(3.16a, b) 


Here the coefficients Aj , A 2 are defined by 

2a 20 XH 2 


\ _ r n , ^ 20^11 2 1/r v sy , «20 , 2X n \H 2 , 

2 1 20 2 + 3(a 20 X„A + /W) I/[ ” G2+ a| (1 +^ 0 + 3(a 2 „X„A + A„^)l’ 

(3.17a, b) 

where for n = 1 , 2, G n ,H n are defined by 


G n = 


5A 


a noX n + /? n o 3(a n 0 AT n A + /3 n op) 1 


Al'(Zon) Xo(M' 


with 


6 n = 


— i sfi 


nO 


(a n oX n A + Pnop) 3 


010 = 0 . 
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Finally the constants ai, bi, a 2 , b 2 are defined by 

= 

B V ,\X„ /“ <r, «';,<(( - 2 tPt' u fcc)A' n ’ 

. = 0<« ir r,i' 4 d( 

1 /} I0 AX„f 0 °°<;r 1 w; I d(-2iw n (co)x„’ 



“ 2 A A„ /“ Cr 2 [a 2 „! 7; 2 + - 2 t [a 2 „t / 12 + ’ 

t = / 0 °°r 2 [a 20 ^ +/? 20 $;]< 

2 AA„ /“ Cr 2 [a 2 „c; 2 + &„ w; 2 ]dc - 2i[o 20 C/ 12 + ’ 

where ri and r2 correspond to ‘r’ in ( 3 . 15 ) with (o:,/?,fl) = ( 0 , /?io, flio) and (a,/?,n) = 
(<*20, #20,^20) respectively. 


4. The generalization and solution of the amplitude equations 

Firstly we note that the three-dimensional wave can also be ‘de-tuned’ by varying di 
by an amount ef/? from the neutral value. In that case (3.16) becomes 


-^ = \jA~a x A\A\ 2 -b x A\B\\ 
dX 

- Sr = As )b + A 2 XB - a 2 B | B | 2 -b 2 B | A | 2 , 


(4.1a, b) 


where A3 is defined by an expression similar to (3.17a). Secondly we note that (3.16) 
apply in a C2 neighbourhood of the position where the three-dimensional wave is neutrally 
stable. Following Hall and Smith (1984) it can be shown that (3.16) apply over a longer 
length scale if X 2 X in (4.1b) is replaced by X 2 X where X is then treated as a constant in 
the amplitude equations. This result can be found directly from (4.1) by letting X — ► 00 
and introducing a length scale shorter than X in order to retain the derivative terms. 
The resulting amplitude equations have a ‘quasi-parallel’ nature and correspond to the 
calculation of Smith (1979b). 

We now define p and <7 by 


p =| A | 2 , cr=| B 


in which case (3.16) and the generalization of this system for X » 1 can be written 


P\ — 2p{Ai r /3 ®lrp ^lr^}, 

<7y -— 2(7 { A2r-X^ ®2 ^2rP), 


(4.2a, b) 
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(4.3a, b ) 


and 

Px = 2p{Ai r ^ &lrp b\ r <7 }, 

(Jx ~ 2<t{A 3r /? + A2r-? — a 2r fr — b 2r p}. 

The precise nature of the solutions of (4.2), (4.3) depends sensitively on the constants 
appearing in these equations. We shall see in the next section that the constants a\ r and 
a 2r are positive almost everywhere so we first discuss such a situation in detail. In fact 
ai r is always positive and this result is entirely consistent with the finite Reynolds number 
calculations of Hall and Malik (1986) for the two-dimensional mode. 

A matter of some importance is the question of whether p or a in (4.2) or (4.3) can 
become infinite at a finite value of X. This would mean that three-dimensionality could 
destroy the stable equilibrium states of Hall and Malik (1986). We seek a singularity of 
either system as X — > Xo by writing 


P = 


P o 


a = 


(X 0 -X) 

fro 

(X 0 -X) 


+ 

+ 


in which case po , fro satisfy 


2 — a lrP0 bi r <7() , 
- = —a 2r o o — b 2r po, 


(4.4a, b) 


and po and <7o must of course both be positive. It follows immediately that no such 
singularity is possible if a,\ r . a 2r , &ir and b 2r are all positive. In fact it is easily shown that 
with ai r and a 2r positive the only case when the singularity can occur is when 6 lr and b 2r 
are negative and 

ai r a2r ^ birb 2r - (4.6) 


This condition effectively identifies an important class of three-dimensional waves which 
can have a significant effect on the two-dimensional equilibrium states of Hall and Malik 
(1986). In order to see why this is the case it is necessary for us to discuss the solutions of 
(4.2) and (4.3) in more detail. We continue to discuss the solution for the case when ai r 
and a 2r are both positive. 
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In fact we begin with a discussion of (4.3) and return to (4.2) later. It is easily shown 
from (4.3) that p and a have the possible equilibrium states: 

a. p = cr = 0, 

b. p = Air/ja^ 1 , a = 0, 

C. p = 0, O — [\ 3 J + \2 r X]a 2 r , 

d. P = [A i r /? — ftirCrJaj,. 1 , <7 = {Ai r /?&2r ~ [A 31 "/? + ^ 2 rX]a\ r } / {bi r b 2r — ai r a 2r } . 

(4.6a, b, c, d) 

The solutions b and c correspond to ‘pure’ 2D and 3D modes respectively whilst d is a 
mixed mode. If the detuning parameters /3 and (3 are held fixed whilst X is varied we can 
determine the evolution of the equilibrium amplitudes as the disturbance develops away 
from the attachment line X n = 0. Note that A 2 r from (3.17b) is negative. The stability of 
the different equilibrium solutions can be checked by a routine stability analysis. Before 
discussing the nature of the solutions we note that in all the cases we computed, a i r and 
02 r are almost always positive in which case nonlinear effects are stabilizing if either the 2D 
or 3D mode exist separately. We further assume that the detuning parameter (3 has been 
chosen such that \\ r j3a\ r > 0 so that in the absence of a 3D wave a stable finite amplitude 
wave propagating along the attachment line is possible. If a\ r and a 2r are positive then 
there are four possible combinations of signs for bi r and b 2r . The bifurcation properties 
for these four cases are summarized below: 

CaSG I dir, 02r ) b\r, b 2 r ^ 6* 

The different possible solutions in this case axe shown in Figures 3a, b for the ‘sub- 
cases’ b\ r b 2r > a\r a 2 r and b\ r b 2r < ai r a 2 r respectively. Sufficiently fax upstream we see 
that only the pure 3D mode is a possible stable mode whilst sufficiently far downstream 
only the 2D mode is a possible equilibrium flow. In the case b\ r b 2r < a\ r a 2r there is a 
short interval where the mixed mode is the only possible stable state. 

Case II ai r ,a 2 r > 0, b\ r ,b 2r < 0. 

The solutions in this case axe shown in Figures 3c, d for the ‘sub-cases’ b\ r b 2r > a\ r a 2r 
and bi r b 2r < air« 2 r respectively. In the first case the only stable solution is the 2D mode 
beyond the position where the mixed mode bifurcates from it. However, a phase plane 
analysis shows that a sufficiently large disturbance to this state is unstable. Thus there 
is a threshold type of response where a small disturbance to the 2D mode decays whilst 
a sufficiently large one will grow. The size of the ‘sufficiently large disturbance’ decreases 
to zero as X decreases to the point where the mixed mode bifurcates. Before this point 
there are no stable modes and any disturbance will grow, in this case and the threshold 
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amplitude case the growing disturbances terminate in the finite X singularity discussed 
previously. 

Case III ai r ,a 2 r,hr > 0, h r < 0. 

Here the situation is as illustrated in Figure 3e. Dependent on the value of X either 
the mixed or 2D mode is stable. A phase plane analysis shows that each stable state 
is stable to an arbitrarily large disturbance so there is no threshold amplitude type of 
response. 

Case IV j CX 2 7* 5 0, &2r ^ 0. 

The situation is now virtually the same as Case III except that the mixed mode loses 
stability to the 3D mode when X decreases so that the 3D mode is stable as X — ► — oo. 
Again there is no threshold amplitude type of response at any value of X . This is illustrated 
in Figure 3f. 

Thus we see that apart from the case ai r ,a 2r > 0,6i r ,&2r < 0 with &i r f> 2 r > air« 2 r 
there is always a stable equilibrium state available at any value of X. Furthermore in the 
latter situation at sufficiently negative values of X the stable state is never the 2D mode. 
However, as the disturbance develops with increasing X ultimately only the 2D mode is 
stable. In the exceptional case a sufficiently large initial disturbance will terminate in a 
singularity at a finite value of X. 

We now turn to the case where a ir and a 2r are not both positive. We shall see in 
the next section that this situation is unusual and occurs when the constant a 2r becomes 
negative so that nonlinear effects destabilize the three-dimensional mode. The situation 
in this case can be investigated following the previous discussion. The main result is that 
(4.1) then always permits a solution which becomes infinite at a finite value of X. The 
singularity has the same structure as that discussed above with the only change being 
that, dependent on the other constants, it is possible for B alone to become infinite. The 
equilibrium solutions of the amplitude equations and thus instability characteristics can 
similarly be investigated for the case a 2r < 0. Here the three-dimensional mode bifurcates 
to the right and is always unstable. In some situations the mixed mode exists and it 
is possible for the two-dimensional mode to be stable to small perturbations. However 
sufficiently large perturbations always destabilize the flow so that we conclude that when 
a 2r < 0 the presence of sufficiently large amplitude perturbations will always lead to the 
finite-X singularity being set up. We conclude that there are just two situations where the 
ultimate state set up after a wave interaction between two and three-dimensional modes 
will not be a stable two-dimensional mode. These exceptional circumstances correspond 
to when ai r ,a 2r > 0,f>i r ,& 2r < 0 with ai r a 2r < &ir& 2 r or whenever a 2r < 0. 
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A similar type of discussion for ( 4 . 2 ) is not possible because there are no equilibrium 
states for this system for all X. However, for large values of X it is easy to show that 
there is a solution with p = Ai r pa^, a = 0 and that this solution is stable. There 
are no other equilibrium states so that, unless limit cycle solutions of ( 4 . 2 ) exist, or a 
singularity develops we expect any initial disturbance to evolve into a pure 2D mode at 
large X. Numerical investigation of ( 4 . 2 ) showed no evidence of limit cycle behaviour and 
that in the exceptional case a finite X singularity develops and the 2D equilibrium state is 
then never set up. It remains for us to discuss the values of a\ r , a,2r, hr, hr found in our 
calculations so that the above results can be applied to the instability of attachment line 
flow. 


5. Results and discussion 

We have seen in the previous section that the nature of the solutions of the amplitude 
equations depends crucially on the constants ai r , &i r , a2r, &2r- These constants can be 
found only after the differential systems for the fundamentals, adjoint, first harmonic, 
mean flow correction functions have been solved numerically. These systems were solved 
using finite differences in the manner described in Hall and Smith ( 1984 ); the reader is 
referred to that paper for a more detailed description of the method. It was found to be 
convenient to map the region 0 < ( < oo into [0, 1] using the transformation 

2 _i , 

77 = - tan C, 

7T 

which aids the convergence of the velocity field at large (. The other significant difference 
between our calculations and those of Hall and Smith is that here the spanwise momentum 
equation has a solution with the velocity component tending to a constant rather than 
decaying algebraically to zero. In order to illustrate how this can be taken into account we 
consider the equations for (U21, V21, W21, P21). By combining the X* and Z momentum 
equations we can show that F = 020^21 + and G = pJJ2i — ^X n W2\ satisfy 


F"' - [-ift. 20 + i(\X n a 20 + nPio )(}F' = 0, 

G" — [— if22o + i(\X n a 2 o + Vp2o)(]G = *(<*20^ — A-X n /?2o)-f2i • 


The first of these equations is to be solved such that F( 0 ) = F'( 0 ) and F — ► constant 
when £ — * 00 whilst the second equation is solved subject to Cr(0) = 0, C?(oo) ~ £ -1 . Thus 
the combination p,U2i — A X n W2i decays algebraically when £ — > 00. Having solved for F 
and G they can be combined to determine U21 and W21, then the equation of continuity is 
solved to determine V21 . The equation for the first harmonic functions can be integrated 
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using the same procedure. Finally in our discussion of the numerical scheme we note that 
the convergence of our scheme was checked when appropriate by varying the step length 
over Ccc the approximation to oo in the £ direction. 

The constants <*1,61,02 an d 62 were calculated for X n = 0.1, 1., 5. and 10. The results 
are normalized by making a\oU[\ + (3ioW[\ and <*20^12 + #20 W(' 2 both equal to unity at 
£ = 0. Some typical value of these constants are shown in Table 1 below. We see that 
it is possible for either of the two exceptional cases of the previous section to occur. In 
Figures 4 and 5 we have plotted the neutral values of the a 20,^20 and indicated where 
the exceptional cases occur. The first exceptional case with ai r ,a2r > 0 is denoted by • ■ • 
whilst the other exceptional case is denoted by — — . 

We see that at X n =0.1 an interaction of the two-dimensional mode with the three- 
dimensional mode with a//3 >~ 1.27 will cause a singularity in the disturbance amplitudes 
to occur. Thus at X n = 0.1 three-dimensional waves propagating at an angle of more 
than about 50° to the attachment line will cause the catastrophic breakdown of the two- 
dimensional mode. 

A further band of modes with <*20 < 0 which leads to the first exceptional case is also 
seen to exist. These correspond to low frequency three-dimensional modes. In the limit 
as «20 — ► 0 these modes have zero effective shear stress and correspond to the stationary 
viscous crossflow modes of Hall (1986) and MacKerrell (1987). We conclude that near 
the attachment line the stimulation of oblique waves propagating at an angle greater than 
about 50° or the stimulation of the viscous crossflow modes of Hall and MacKerrell will 
cause a new larger amplitude disturbance flow structure to develop. 

When X n = 1 only the destabilizing band of wavenumbers corresponding to the low 
frequency modes remains and the interval over which they exist has decreased. However 
when X n = 5 the stationary viscous crossflow modes become subcritically unstable so that 
the stationary viscous crossflow modes cause the finite X singularity to develop at almost 
all of the possible negative values of <120 • In addition there is a very short band of oblique 
modes propagating at an angle of about 80° to the attachment line which leads to the 
singularity being set up. This band of unstable wavelengths no longer occurs at X n = 10. 
but the stationary viscous crossflow modes are now subcritically unstable for almost all of 
the possible values of a 2 o with a 2 o < 0. 

Without prohibitively expensive numerical calculations we cannot confirm that the 
results discussed above show the overall trend of the possible interactions when X n in- 
creases. In fact some further investigation showed that the small band of destabilizing 
oblique modes at X n appears and disappears as X n varies. However our calculations do 
suggest that at small values of X n there is a wide range of possible oblique modes and a 
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small band of low frequency modes which, if excited, will cause a catastrophic breakdown 
of the disturbance flowfield. Further away from the attachment line the oblique modes 
become less important and it is the lower frequency modes which become the dominant 
mechanism. 

Clearly our analysis cannot predict what kind of flow will be set up once the singularity 
appears. However we note that other modes, notably the inviscid stationary crossflow 
vortex mode of Gregory, Stuart and Walker (1955) might then become important. 
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The neutral eigenvalues satisfying (2.24) 

The neutral eigenvalues satisfying (2.24) 

The imaginary part of ctg satisfying (2.24), as a 
function of X for a range of values of 0 . 


The equilibrium solutions for a lr> a 2 , b lr< b 2r > 0 
with b lr b 2r > a lr a 2r . 
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